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I. Prefaoe 

The mathematical equations in this document form the basis of a computer program 
used to perform an analysis of the errors that oan be expected in a set of post- 
flight aerodynamio parameters. These errors are due to inaoouraoiea in the Shut- 
tle best estimate trajeetory (BBT) and in the meteorological data to be provided 
by the national Heather Service in support of the OPT flights. 



XX. Introduction 

The t « to »> of this nfort to to dm bon to start with o give* 
a tat* vector, X, and its associated error covariance matrix, C^, 

sad free these calculate th* parameter vector, I, and its associated 
error covariance matrix. 

He have 



shuttle (vehicle) BCI position vector 
(meters) * (a) 



shuttle BCI velocity vector (a/ sec) 


IMU platform aisallnsaent angles (radians) 
■ (rad) 


—ATM, TOP 
P 


wind velocity in topodetic coordinates 

(«/•«) 

" " 3 

atmospheric density (kg/m ) 


T 


atmospheric temperature (*K), 


and C x is a 14x14 matrix. The source of X and will be the 
following: 


X 


I — 

M 

id 

(*! 

i 


^ECI 




• 

^ATM.TOP 


P 


T 



postflight best estimate trajectory (BET) 


National Meteorological Center (NMC) 


the upper left hand 9x9 submatrix of by the BET, the lower right 
hand 5x5 submatrix of C x by NMC, and the rest are zeros; l.c., 


2 


... . 






( 


\ 




{ 

{ 




The parameter matrix, Z_, la defined aa 


Z 



vehicle airapeed (m/aec) 

fllghtpath angle vrt etmoaphere (rad) 

aslmuth angle vrt etmoaphere (red) 

slfltude above the Flacher ellipsoid (m) 

2 2 

dynamic pressure (newtons /a ) ■ (nt/m ) 
equivalent velocity (read by pilot) (m/aec) 

Mach number (-) 


V» 

CO 

a 

s 

0 ' 

q« 

q0 

? 0 ' 


hypersonic viscous parameter (-) 
angle of attack (rad) 
sideslip angle (stability axis) (rad) 
slldalip angle (body axis) (rad) 
pitch dynamic pressure (nt rad/m ) 

2 

yaw dynamic pressure (stability axis) (nt rad/m ) 

2 

yaw dynamic pressure (body axis) (nt rad/m ) 


All calculations In this report will be In the metric system. 
Ill . Calculations 

C z Is calculated by using the Identity 



3 





37. 

vhar* P ■ -r la the 14x14 matrix of partial* 

0A 

(C, la also 14x14). So our problem la raducad to calculating l 
and P. 

notation : EC1 • Earth Centered Inertial 

EF * Earth Fixed (non- Inertial) 

TOP - Topodetlc (non-lnertlal) 

V - Vehicle 
ATH • Atmosphere 
B - Body 
P ■ Platform 
M? - Misallned Platform 
Coordinate Transformation : 



*EF - < T4 > —EC! 
where (T4) - (ROT) (RNP) , 


(ROT) 


coa (Dg (T-T q ) 
-sin u B (T-T q ) 
0 


a In «o E (T-T o ) 

cos u_ (T-T ) 
E 0 

0 


0 

0 

1 


( 2 ) 


0) 

( 4 ) 

(5) 




j 

i 

i 




(accounts for rotation of earth) , 

(RNP) is the rotation, nutation, and precession matrix, T q is 
its time tag, and 


w E ■ 7.2921159 x 10 -5 ( “ radians the earth turns per sec) 

Since the topodetic coordinate system is centered at the vehicle rather 
than the renter of the earth it is necessary, when transforming points, 
to first translate and then rotate. However, for vectors where only 
direction matters (e.g., velocity, acceleration, etc.) It is only 
necessary to rotate coordinates to get topodetlc vectors to EF vectors 
and back. Therefore, the following transformation is used for such 
vectors: 



) 


4 



where 


^POP 

■ <«)»„ 



(6) 

ft 


-ala 4 coa X 

-aln 4 aln X 

coa 4 


(T5) 

m 

-aln X 

coa X 

0 

(7) 



-coa 4 coa X 

-coa 4 aln X 

-aln 4 



and 4 (geodetic latitude) end X (geodetic longitude) are calculated 
aa follower 


aet e * 1/298.3 (elllptlclty of earth) 

Rg ■ 6378166 n (radlua of earth) 

B - .0067 (Initial B) 
o 

and let 

e(2-a)Rg 

V°4f + Y CT )/<B i +1)2+ (1 - )2z ep 

for 1*0, 1, 2, 3 

define 


B " V 


Then 


X - arctan ! 


4 ■ arctan 


(ft) 


Zgp(Bfl) 




EF 


( 8 ) 

( 9 ) 

( 10 ) 


(ID 


( 12 ) 


(13) 


( 14 ) 
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Now Co calculate Cha actual velocity of tha vahicla over tha ground 
va work In EF coordinates • 


( ) 


5 ef “ (T4) hn 


ao 


O) 


*EF " < T4 > ^CI + < T4 > *ECI 


and expressing this In TOP coordinates 


Stop h (t5) ^ef 


(15) 


(16) 


or 


*TOP - (T5) (T4) R SC1 + (T5) (T4) ^ 

or defining (see equation 20 for (T4)) 

(T6) - (T5) (T4) 

(17) - (T5) (T4) 


(17) 


(18) 


ve have 


Stop “ < T6 > hci + (T7 > hci 

- £(T6) (T7) (0) (0) 


3xl 

matrix 


3x3 

matrices 


(0)] X . 

1 t 

3*2 14xl 

matrix matrix 


(19) 
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In the above 


(T4) • (EOT) (EH?) 


-sin Mg (T-T e ) 
■ « s -coa « B (T-t 0 ) 


coa Mg (T-T o ) 0 

-a in «g (T-T 0 > 0 


So froat the above ve have 


Ry | ATM TOp s vehicle TOP velocity vrt ataoaphera 

A 

a a 

“ ^rop “ ^Th.TOp 

- |CT6) (T7) (0) (-1) (0)1 X 

1 ' t J t 

3x1 3x3 3x2 14xl 






m 

[ (T6) 

<T7) (0) 

(-« 

(0) ] 

♦ 1 

. Tx (Ti> 


>> 

*x ( - 7) 


3*'** A* 


3x2x14 


• £ (T6) (T7) (0) (-1) (0) ] 

+ [( ■Sc (M) ) < w> ( Fx < t5) ) (T4> W (0> <°>j * 

- [ (T6) (T7) (0) (-1) (0) ] 

+ £(^(T5))<T4) (~(T5»(T4) (0) (0) (0) j X 


3*3x14 


3*3, *14 3*3. *14 3*2, *14 14; *1 

1 a I 1 1 

3*3j*14 


- [(T6> (T7) (0) (-1) (0)] 

+ (^-<T5))(T4) (^<T5))(?4) (0) (0) (0) 

-V _ _ _ 4_ . L 


^(•^(T5))(T4) (-^(T5))(T4) (0> (0) (0)j X || 


8 ^v|atm.top 

3X 


(T6) (T7) (0) (-1) (0) 








3x3, 


3x3. 


3*1. 


3x3, 


3x3, 


vhar* (T6) and (T7) ar* defined in aquation (18), 
(T4) in aquation (4), 

(T4) In aquation (20), 

(T5) in aquation (7), 

X In aquation (13), 

4 In aquation (14), 


£(T5> - 


a in 4 aln X 
-eoa X 
coa 4 aln X 


-aln 4 eoa X 
-aln X 

-eoa 4 eoa X 


(0) 


IX 

lx 

A 

3*®i 

T 

l V l 2 

f 

1 2 *14 




3x14 



(0) 

j?h 

li 

IX 

3«8 l 

T 

u i M h 

f 

1 2 *14 


0 

0 

0 


( 22 ) 


(23) 


ir 


(T5) - 


-eoa 4 coa X 
0 

aln 4 coa X 


-eoa 4 *in X 
0 


-aln 4 
0 


aln 4 aln X -coa 4 


(24) 




ul_ . 

>h. 


$ ■&»» - m 0 


kfr™ - + r 2 + T 2 ) 1/Z 

pT^ **"*• " «» 

r-i^ttH) 2 


2,1/2 »■_ 
Bt 1 »T_ 


L . c f ot „ 2 ♦ T^V'Vl) * <I „ 2 * V >1/! ft] 


V + V + *bf 2(b+1)2 


44... !zSL 

2*0, « 


- If- • H- • H- < T4) (0) (0) <0) (0> 

**EF ™St 8Z E t I 


1. V„ - vehicle airspeed - l»v | ATMiT01> j < ' <> 

3V T _ W T ^lATM.TOP 

S - a» . ax 

^.’]ath,top 


. fe I ATM. TO?) a 3y|ATM.TOP 

I^VjATO.TOpI 1" 

t r 

1*14 1*3 3*14 [see equation* (21), (22)] 

aV T 

Thl« 1*14 matrix, la the flret row of tba matrix of partial*, 

P [aee equation (1)7- 


C) 




^~^-^5^^WS«ar i III | > 1^,1 i; y W JJ j^ i P^Wji i' |h_ | i_ ^ l l J ' i ni ' I 'll Il l 'll , I 


m r wy* ww !?5P w , M? '^w#! 



,TOP 

X-QJ. (NORTH) 


Figure 1. 


2. y_ » flightpath angle wrt ataeophere 


■ arcs In 


« 

[ ^1 ATM, TOP 1 

I^V|ATM,TOpl J 


(rad) 


% 

3X 

* 


3r n 


^ I atm. top 


1x14 


3 ^v|atm,top - 

( _1<, I^v|aim,topI S x v|atm,top + Y v|atm,top) 

[ _X v|ATM,TOP ^vj ATM, TOP* " Y V I ATM, TOP ^V| ATM, TOP* 
V] ATM, TOP * | ATM, TOp] 


*2 

X 



3. ^ “ azimuth angle wrt atmosphere 


arctanpl*™.™?] 

- *v Iatm TOP -* 


(rad) 


(29) 


(30) 


(31) 


(32) 
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ax 


3? V AW, 


ISytATW.TOP 

ax 


TOP 


-J > — T2 

X v| ATM, TOP + * V I ATM, TOP 

[-*v|atk,top* *v|atk,top» °] 


ATM. TOP 


1x14 


ftX 

T 

3X14 


( 33 ) 


(34) 


4. H ■ altitude above the 1960 Fischer ellipsoid 




+ Y^)/(W-1) 2 


+ 



<■) 


where e and B are defined above (see page 3) . 


(35) 


For the error calculation we use the following approxiaatlon for H, 
which is accurate to within 1 meter for altitudes up to 200,000 meters! 


B - ISeci! - 


yi-e)lSBci I 

V ( l- e)2( 4p + + 2 


2 

EP 


but we know that jR^jl » |Rgp| . The equation then becomes: 

H “ I®*,! 


R (l-e) 

C 


fT- 


•> 2 <V * V 2 > + z a 
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n«t 

vu 


fl-.airi-<l-e> 2 > * 2 

*SF 

T5T r 

[(l-e) 2 (X 2 KF + T 2 gy) + t 

an 

v 

EF 

(1-e) [1- (1-e) 2 3^ Z 2 


(1-e) [!-(!-«) 2 JS^OT + T 2 ) 


|^(l-e) (X 


and for small e we have 


* WZ T 1 + 2e 


*E EF ^ ^Ep 


and in fact e ■ 1/298.3 Implies that 


,.!« 2 g . JL ad , r e (x ef + r 

2 | R 1 3 149 * and 2e ~T ^3 

'-KP ' I- Kp! 


( X EF + ^EF^ 1 


so th&t as claimed above the second term Is small compared to the first 




Therefore, 

IH _ 3H 8 S» 


i 

3 


3H as 31 
ax » ’ 8Y ef ’ 3Z ef 

1*3 


< 

f 

- 


(tt) (0) (0) (0) (C) 


3*14 


3H 


where — taken from equation (38) is dose enough. 
3 ®BF 

5. q « dynamic pressure ■ yp V* (nt/a ) 


ia -f 1 i£ + 2 3 T 1 

3X " q [p » ^ 3TJ 

- ( (0) (0) (0) (0) 


1*14 


1*12 


2 *Tl 

d.O) ) + arj 

I f 

1*2 1 * 1 * 


8V W 


(see equation (28) for yjr 


Here we neglect terms like 


j£_ 


3 ^eci 


6. V, 


EQ 


3V t 

since they are small compared to , etc 


equivalent velocity ■ (m/sec) 


where p q - atmospheric density at sea level 


- 1.2250 kg/m 3 (1952, standard atmosphere) 


!Is - 5 a 

3X 3X 

2q 

[see equation (42)] 


(41) 


(42) 


(*3) 


(44) 


( 45 ) 
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and 



T ♦ 122.1 x 10~ ( y^ 

T» + 122.1 x 10 -(3/T '^. 


where 

T’ - 726.97 + 



,468T + 3.63921 x 10 

3M m 

00 

1 

3C^ 

3X 


3X~ 

3M - 

1 

r 3c : 

3X 

+ 2C r 
00 

[ 3T 



3T . ^ 3T , 1 

3X W 3X J 


1 Fi_ !!i , i i£ i ieI 

2 j_ V T 3X P 3X " v 3XJ 


3C * 
00 

W 



1.5 

T 


+ 


1 + (log 10) • 122.1 x 10 


“(5/T) 


V* 

I 2 
s w 


T + 122.1 x 1C 


-(5/T) 





1 + (log 10) • 122.1 x 10 


-(5/T') 


..2 


T' + 322.1 x 10 


(5/T’) 


3T* 

3X 


.463 § + 2 


(3.63521 x 10‘ 5 ) V T 



3m 3p 3T 
3X " 3T 3X 
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Mor« Coordinate Trana formation. . 
Recall equation (15)} 




k? * hci 4 CM) 

MW1 


(IS) 


where (T4) an? (T4) are defined In equatione (20) and (4), respectively. 
Alao we have 




TOP 


(61) 


where (T5) 1# defined in equation (7). 
So that 


*»!«*,» ■ «*> ha + <«> ha 


05)1 hm.ra 


(62) 


From here we can go to body coordlnatea: 


Sv|*»,B - <«> m <y (M) 1 i, |Aniilp 


(63) 


*"• (W) T rotate, back to ECI; (I p ) t, a ton.Ua, aatrla that tr^foti. 


from ECI to nominal platform coordlnatea; l.e., 

" C T p ) ^ CI ; 


(64) 


"orSS“£” ST ” 1 PUtt °" “O' 41 "'" “ ■«»!. ala.llnad 


«>v 


hr 

where T Is given by 


(65) 


(T) - 


C ,C , 
p2 p3 


C pl S p3 + S pl S p2 C p3 


S pl S p3 ~ C pl S p2 C p3 


" C p2 S p3 


C pl C p3 " S pl S p2 S p3 


S pl C p3 + C pl S p2 S p3 


P2 


-S ,C _ 
pl p2 


C t C , 
pl p2 


(66) 
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(bar* C pl - eoa » pl , S p2 • ala 6^, an d similarly far the others) 
and a Inca tha angles « pl , 6 p2 , and # p3 are small , we have for arror 
analysis purposes : (T) «■ (T'> where 


(T') 


1 

•n 

-*n 

p3 

1 

•pi 

P2 

‘V 

1 


(67) 


(9p is part of tha given state vector X); and where (Q) 
the nlsallned platform coordinates to body coordinates; 


transforms 
l.e. f 


^v|atm,b “ 5v|atm,mp* 


( 68 ) 


where Q Is a transformation matrix formed by rotation of tba platform 
about Its axes by their respective almbal aneles. 



T * 

* ti 


5v|atm,b " (Q) (T) <T p ) <T4> ^|atm,ep 


(Q) (T) (T p ) (T4) J 


[<«> bci + (T4) ha " <t5)T ^MM.TOp] 


<V [' 


«> (!) (T„> | («) T <T4) ig,,, + V, - («) W i*a, TO p] 


or defining 


(69) 


(T8) - (T4) T (T4) 


(T9) - (T4) T (T5) T 


(70) 


20 


w hm 


«J) (T) (T o > [ <T8) ^ - (T9) | 1 

“ w MiAfiyr j 


*v|aw,b 

m 4 Mln * h#rt *' la*t«d of T for tbo orror enely.lo: 

■ W) (T') CT p ) £<T,) ;i) (0) (_IJ) (0)1 

» « ( 0 ) ] 

- %*■* la givaa byi 


09 

**ECI 




(71) 


(72) 


’ W) (^) ‘V [<V *« + - <”> »,».«»] «» 


whore 


ill 

36 


3X» 

W~ 


3T* 

dTT 

P3 


r 


0 

0 

0 

0 

0 

1 

" 0 
-1 
0 


0 

0 

-1 

0 

0 

0 

1 

0 

0 


0 

1 

0 

-l" 

0 

0 

0 

0 

0 


( 74 ) 
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10. 0 * sideslip M|li (stability uli) • arctan 




*i *2 • 

ag . SIatm.b + S Iatm.b t tIaim.b 

V^lAIK.. * iJlUM 


( 




"Si AIM.*, 


*2 «2 • • # *2 *2 """ 
*V jATM,B + SjATM.B S|ATM,B SIaTX.B + S|a» 


■ aln 0 cot 6 


( ;S|ATM.B 1 ‘SjATH.B 

VS I ATM, B + S|ATM,B T V|a1M,B *V|aTM,B * 7 


li _ 30 8 2 vIatm.b 

9 I as « 

35 v|atm,b 


V xi vU tm.b + zZ v|atm,b ( x %|atm,b • t2 v1atm.b + x %|atm,b) 

^ X V|ATM,B S|atM,B* X J|ATM,B * 2 V|ATM,B* “*V|ATK,B Z V|ATK,i 

^lATM.B 


IX 


11. 8' “ aide* lip angle (body axis) 

/s1atm,b\ 


arctan 




< 33 ) 

( 84 ) 

( 85 ) 

( 86 ) 

( 87 ) 

( 88 ) 




<T*> « (T5> (14) 
(T7) - (T5) <T4) 


( 18 ) 


*^r|ATH.TOF 

a 

T 

3*1* 


- j(«) (T7) (0) (- 


3*2 


r i 

[lx- («)] (T4)j ^ (T5)] (T4)j (0) (0) (0) X 


Si 

3T 


3*14 


t J t 

14*1 1*14 


[[£ <«>] <«>] [[i <”>] <»>] 


(0) (0) (0) 


r 


li 

ax 


3*14 


14*1 1*14 


ivi*™,iop - [<”> '”> «» <-» <«>] i 


o 

Y R 

l b 


1.2250 kg/m J 

401.872524 m 2 /*K-Bec 2 
32.765 a 

1.458001 10“ 6 1 l, lt/«ec 

T + llO 2 


(22) 

( 21 ) 

(44) 

(47) 

(50) 


- 

(T p ) - 


REFSMKAT 

3*3 matrix such that 
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T is given after the flight. 
P 



C P2 C P 3 

C pl S p3 + S pl S p2 C p3 

pl p3 pl p2 p3 

(T) - 

“Vp3 

C pl C p3 ~ ®pl S p2 S p3 

S ,C , + C ,S „S , 
pl p3 pl p2 p3 



-S .C , 
Pl p2 

C pl C p2 j 


(C p l = cos 

8 pl* S P1 5 8l “ V* etC<) 



Q - 

(T8) - (T4) T (T4) 

(T9) - (T4) T (T5) T 

V - (y [<»)«„!♦«„,- (19) S» TO>TOP ] 

\ 

0 0 

0 1 

1 0 

0 -1 

0 0 > 

0 0 


5v|atm,b " <« < T > - 


3T’ 


36 


Pi 


0 

0 

0 


3T 1 

36 


pl 


0 

0 

1 



3*3 netrlx 

(see discussion, page 18 
after equation (68)} 
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( 78 ) 

(81) 

(83) 

(85) 

(87) 


fl 

3X 


aw * T * 

. I . jy I atm, top) ^vIatm.top 

- Z 1 3X 


(28) 


^2 % (-^lATM.TOP^ViW 

» ‘ST " p . 


» Y* • v 

,TOP Vj ATM, TOP/ 

r • 

L V|ATM,TOP ^V J ATM, TOP * ~ Y v|ATM, 


X 2 ‘2 I 

VjATM.TOP + Y V | ATM , TOPj *h\kTK. 


TOP V | ATM, TOP* 

TOP 


3X 

7 

3X14 


(31) 


ax 


a*. 


R 


ax ■ 75 


X V| ATM, TOP + Y VjATM.TOP 
[" Y V| ATM, TOP* X V | ATM, TOP * °] 9 -i 


VI ATM. TOP 

ax 

♦ 

3X14 


(34) 




31 






( 84 ) 


!!a 

. ifigl m 

><w 

■ i. 

M » 


i 

W 5 

n 

ax 

ax 

3 

ax 


9 

JX 


. AML - 

*<VlO> 

• Z. 

“10 

+ 

l 

«5 

ax 

» 

ax 

5 

ax 


10 

ax 

3Z 14 

- AMU . 

•w 

• V 





ax 

ax 

ax 

*5 

ax 

▼ 

Z 11 

ax 


( 86 ) 


( 88 ) 


MOTSs Th« IMS program dou not utilise tha covariance matrix . 

Instead it Inputs the error values for the state vector, 
AX. (1-1, . . .14) and calculates KS8 output errors by: 
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